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We explore a model for a fermionic dark matter particle family which decouples from the rest 
of the particles when at least all standard model particles are in equilibrium. We calculate the 
allowed ranges for mass and chemical potential to be compatible with big bang nucleosynthesis 
(BBN) calculations and WMAP-data for a flat universe with dark energy {0.% = 0.72, Q.°M = 0.27, 
h = 0.7). Futhermore we estimate the free streaming length for fermions and antifermions to allow 
comparison to large scale structure data (LSS). We find that for dark matter decoupling when 
all standard model particles are present even the least restrictive combined BBN calculation and 
WMAP results allow us to constrain the initial dark matter chemical potential to a highest value of 
6.3 times the dark matter temperature. In this case the resulting mass range is at most 1.8 eV < 
m < 53 eV, where the upper bound scales linearly with glffiToec)- From LSS we find that similar 
to ordinary warm dark matter models the particle mass has to be larger than ~ 500 eV (meaning 
glf fiToec) > 10^) to be compatible with observations of the Ly a forest at high redshift, but still 
the dark matter chemical potential over temperature ratio can exceed unity. 



I. INTRODUCTION 

In recent years observations and the high precision 
analysis of the cosmic microwave background radiation 
[l|, of large scale structure [Ih of the abundancies of 
hght elements in the universe [3| and of the evolution of 
the Hubble parameter by observation of the luminosity- 
redshift relation of supernovae of type la [3| , @| have been 
strongly indicating that the universe is flat and mainly 
consists of two unknown components, dark matter and 
dark energy. Dark matter appears to be consisting of 
non-relativistic particles that only interact gravitation- 
ally and perhaps by weak interactions. In addition there 
has been a recent independent support for the dark mat- 
ter paradigm by combined x-ray and weak lensing ob- 
servation of the aftermath of a collision of two galaxies 
clusters 0. 

The nature of the second unknown component, dark 
energy, is even less clear, it is mostly characterized by 
having a negative pressure that allows to explain the ob- 
served late time exponential expansion of the universe. 
In the simplest models dark energy is just a constant 
vacuum energy. Combined with a sizeble amount of cold 
dark matter this constitues the so called ACDM or con- 
cordance model. 

For an alternative to dark energy using a large scale 
local inhomogeneity see e.g. 0, [!]■ Some authors have 
even tried to unify dark matter and dark energy within 
so called Chaplygin gas models d, [13, [HI, [i2| • This ap- 
proach faces the difficulty of having a speed of sound 
increasing at late times which makes it probably unable 
to reproduce the matter power spectrum of density fluc- 
tuations on smaller scales. 

Also there has been the phenomenological MOND ap- 
proach by Milgrom [13] using a modified Newtonian force 
law and the more recent TeVeS theory by Bekenstein [l3| 
that tried to explain the observations without the need 
for dark matter. Still most of the work using these the- 



ories has been done on galactic rotation curves, stellar 
velocities in dwarf galaxies and the Pioneer anomaly. It 
is unclear whether they even succeed in explaining the 
observations at these scales so they still have to prove 
being an viable cosmological alternative to dark matter. 
For example Ref. \Td\ find that a simple MONDian model 
with a low matter density gives a rather bad fit to the 
power spectrum of the CMBR anisotropics when includ- 
ing the third acoustic peak extracted from BOOMERanG 
and WMAP first year data. A quite interesting recent ap- 
proach to explain at least galactic rotation curves without 
the need for exotic dark matter or a new theory of grav- 
itation has been done by Cooperstock et al. [l^ using a 
axially symmetric metric and linearized general relativity 
for disk galaxies. Although the mentioned different ap- 
proaches may be able to remove the need for dark matter 
on and below galactic scales it is highly unclear whether 
they can do so for the CMB, large scale structure or weak 
lensing observations. All in all there is at present no con- 
vincing alternative to the dark matter paradigm that can 
explain all the mentioned independent observations. 

The ACDM model is successful in explaining most of 
the observations although it currently predicts an exces- 
sive amount of low mass dark matter halos orbiting larger 
galaxies in both linear structure formation approaches 
as well as in numerical n-body simulations, as found by 
Ref. [13, [Hi- Another problem of the CDM model is the 
prediction of cuspy inner density profiles of dark matter 
halos in dwarf and spiral galaxies, that was not confirmed 
by observations (see e.g. Ref. [H, [13, IMI ) ■ Instead it was 
found that the core density of dark matter must to be 
rather constant. It is unclear whether these problems can 
be resolved in more realistic simulations/models that in- 
clude for example more effects from the baryonic matter, 
or if these actually show that the dark matter properties 
are different from the standard CDM case. For example 
Ref. proposed self-interacting dark matter as a solu- 
tion to the cuspy-core problem, but Ref. [23| found that 
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the interesting parameter range in which self-interacting 
dark matter could flatten halo cores can probably be ex- 
cluded by other observations. Another solution could be 
warm dark matter (WDM) as studied by Ref . f23| , where 
the core densities are reduced by the higher velocity dis- 
persion of WDM. 

It is usually assumed that the dark matter constituent, 
whether it is bosonic (a candidate would be the axion) or 
fermionic (for example the neutralino, the gravitino or a 
massive sterile neutrino), initially has a negligible chem- 
ical potential (or to be more exact that the condition 
T ^ II — m is fuUfilled) . The reason for this assumption 
is the expectation that the matter-antimatter asymmetry 
in the universe is of the same magnitude for all particle 
species. So spin statistics are mostly ignored for dark 
matter and the usual approach for cold dark matter is 
then a free Boltzmann gas. In contrast to that we will 
use a free Fermi gas in our approach to find constraints on 
the allowed matter-antimatter-asymmetry through the 
related chemical potentials and the resulting mass win- 
dow of a thermal fermionic dark matter relic from cosmo- 
logical observations. We assume an additional conserved 
quantum number beyond those of the standard model to 
keep the particle stable, but this is common in many dark 
matter models. We will also explore possible implications 
for structure formation for such a particle species. 



at the GUT scale and at the time of BBN. The tempera- 
ture of equilibrium particles therefore is shifted upwards 
compared to the temperature of decoupled particles (like 
our dark matter particles) leading to the temperature 
relation: 



T = T, 
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(3) 



The equilibrium relations ji = —p, and T = T will still 
be fuUfiUed as long as both fermions and antifermions 
are ultra-relativistic (where barred quantities refer to an- 
tiparticles), although they are decoupled and no inter- 
actions take place anymore. The transition from ultra- 
relativistic to non-relativistic behaviour of both fluids 
needs to be treated numerically, meaning that we can 
start our numerical calculation at relatively late times 
(i.e. after BBN) and simply add the energy density, 
pressure etc. of the underlying photon, neutrino and 
baryon fluids where their behaviour and magnitude is 
known from standard cosmology, WMAP-data and BBN- 
calculations. Neutrinos are treated as massless for sim- 
plicity. 



II. INITIAL ASSUMPTIONS AND 
MOTIVATION 

We assume a thermal distribution of a fermionic parti- 
cle family with 4 degrees of freedom (Spin i) was created 
asymmetrically (i.e. a surplus of particles over antiparti- 
cles or the other way around) at an early stage somewhere 
around the GUT-scale (approx. 10^^ GeV). The assym- 
metric process then became ineffective and the surplus 
of fermions was conserved leading to a non-zero chem- 
ical potential of these particles. We then assume ther- 
mal decoupling at a later stage while it was still ultra- 
relativistic. Thus the particles will be unaffected by the 
following numerous phase transitions and the subsequent 
reduction of the number of equilibrium particle species 
and their degrees of freedom gi. Assuming overall en- 
tropy conservation from the dark matter decoupling scale 
up to today demands: 
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where 
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i— bosons 



i— fermions 



which allows to write the entropy using just one overall 
temperature (usually the photon temperature) . In eq. ([2]) 
it is assumed that all the standard model particles (or 
whatever other particles may be present at dark matter 
decoupling) have negligible chemical potentials at least 



III. ANALYTIC LIMITS FOR ENERGY 
DENSITY AND EQUATION OF STATE 

The energy density and equation of state (EoS) of a 
non-interacting Fermi gas with non-zero chemical poten- 
tial can be calculated analytically in the ultra-relativistic 
as well as in the non-relativistic limit, an outline of the 
derivations will be given below, a more detailed descrip- 
tion can be found in [2^]. 

The energy density e+ ~ e+e and pressure = p+p and 
net number density n~ — n— n of a free non-interacting 
Fermi gas in chemical equilibrium are given in natural 
units (c = fi = fcs = 1) by the following integrals over 
the momentum k of the particle: 
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exp(^) + l 



exp(^) + 1 



(6) 



Where the energy-momentum-dependence reads: 

E{k) = + m^. Here g denotes the degeneracy 

factor (i.e. in our case g = 2 for spin), not counting par- 
ticles/antiparticles since this has already been accounted 
for by the two separate integrals. 



A. Ultra-relativistic limit 

In the ultra-relativistic limit E k, k the integrals are 
solvable giving us the equation of state p — e/3: 
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by using the subtitutions 
k — fi 



Xi 



T 



X2 = 



for the first and the second integral, respectively. By 
renaming xi and X2 to x and rearranging one ends up 
with the following three integrals (assuming fi = —jl to 
be true as said before): 
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So finally we arrive at the ultrarelativistic limit of the 
energy density of a free Fermi gas with non-zero chemical 
potentials : 



120 



9_ 

87r2 



(8) 



This formula can also be derived using the so-called 
Sommerfeld-expansion for eq. ([?]) , which is a power series 
in {T I , as for m = all terms above second order will 
vanish and one obtains cq. ([5]). 



cases to distinguish. Both limits refer to particles as well 
as to antiparticles, since both fluids are still present in 
the non-relativistic regime. 



1. Non-relativistic degenerate limit 

In the degenerate case with {fi~m) 3> T, only particles 
with a momentum higher or equal the Fermi momentum 
kp can be added to the system since all momentum states 
with lower energy have already been filled with particles. 
Eq. (dl simplifies then to: 
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Using eq. Q one sees that the pressure scales as 



9 f k^dk ^ = oc e5/3 (10) 
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which is a polytropic equation of state of the form 
p^Be-r with.4andi.^(^^)''' (11) 

2. Non-relativistic non-degenerate limit 

In the non-degenerate case (/x — m) <g; T, the Fermi- 
Dirac-distribution function can be well approximated by 
the Maxwell-Boltzmann-limit so that: 



9 f ( \ 

e ~ — TT / dkk^ m exp ) 

27r2 / ^ I 2mT J 



X exp [ ^ ) = mn (12) 
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Using the substitution x ~ „ one arrives at 

^ 2ml 

e « 5 m I — 1 exp I ^ j = mn (13) 
whereas a similar calculation for the pressure gives 

3/2 



/ mT ' I fi — m 



P^9Ty — j exp^ ^ 



Tn (14) 



B. Non-relativistic limit 

In the non-relativistic limit E Ki m and consequently 
e ~ m n (and e ~ mn). But there are still two limiting 



Again, the equation of state is a polytrope of the form 
p.Ce^ with,. ^ and C=(-^)''' (15) 
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Interestingly, this leads to almost the same equation of 
state because the exponential factor D = exp (^^^r^) 
turns out to become approximatly constant in the non- 



relativistic limit since fi — m (x a and T 
we will outline below in more detail. 



which 



Normalizing to oq = 1 as before we arrive at the consis- 
tency condition Upoie > a-Rei in the form 
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(23) 



C. Generalized Chaplygin gas 

A polytrope resembles a generalized Chaplygin gas 
which has an equation of state of the form 
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A 
7r 



(16) 



where A and T are normaly taken to be positive. Then 
the Chaplygin gas can be used as an unified model for 
dark matter and dark energy, see e.g. in [§, [l^, fill, fl^. 
The scaling behaviour of the generalized Chaplygin gas 
can be found by solving 
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which can be deduced from energy-momentum conserva- 
tion. The general solution is : 
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(18) 



where eg is the energy density at a = 1. We will now 
examine the non-relativistic Fermi gas solution in com- 
parison. For r = —5/3 the solution is 
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For negative A the energy density runs into a pole at the 

scale parameter apoie = '\J —A/{eQ '^^^ — A). The energy 
density of a free Fermi gas should not become infinite in a 
finite volume. The way out of this problem is to consider 
that the gas will become relativistic, i.e. the EoS should 
change to that of an ultrarelativistic gas p = e/3, which 
can be estimated in the following way. First we formulate 
the constraint in terms of the scale parameter: can af/ei 
be smaller than apoie'^ The point where both limits of the 
non-relativistic and the ultra-relativistic EoS approach 
each other 
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leads to 
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(20) 
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for the energy density of the transition from the non- 
relativistic to the relativistic EoS. Now we can use the 
non-relativistic scaling to eliminate eRei 



a-Rel 

ao 



(22) 



which can be rewritten to 



(24) 



This condition is never fuUfiUed since A is negative, 
meaning that the gas will become relativistic before it 
can run into the pole. 



IV. CHEMICAL POTENTIAL AND 
TEMPERATURE 

A. Ultra-relativistic behaviour 

After chemical decoupling of the dark matter par- 
ticle/antiparticle distribution from the standard model 
particle background, the net number, i.e. the surplus of 
particles over antiparticles of the fermion distribution in 
a comoving volume, will be frozen in (as well as the num- 
ber of particles and antiparticles separately). This can 
be expressed by the conservation equation 



2,Hn- = 



(25) 



where H is the Hubble parameter, leading to a simple 
dependence on the scale parameter a 
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(26) 



The ultrarelativistic limit of the net number density can 
be calculated as a derivative of the pressure or the energy 
density, using eq. ([8]): 



dp_ ^ lde_ j,2 _9_ 3 



(27) 



showing that the net number density goes to zero for a 
vanishing chemical potential as expected. Furthermore 
since the net number density has to scale like for 
any value of T and fj, we can deduce the ultrarelativistic 
scaling of the temperature and the chemical potential: 



T (X a 



so that 



and 



and 



p (X a 



(28) 



(29) 



like for standard relativistic particles (with ^ = 0). This 
is of cause only true as long as the chemical potentials 
and temperatures of fermions and antifermions are well 
above the mass of the particle. For a Fermi gas in chem- 
ical equilibrium with a surplus of particles that becomes 
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non-relativistic, particles and antiparticles annihilate un- 
til all antiparticles are gone (exponentially suppressed 
with respect to particles, since p. —>■ —m) while conserv- 
ing the net number per comoving volume a? ■ n~ . In our 
case not only the net number, but also the numbers of 
particles and antiparticles are conserved separately, forc- 
ing both chemical potentials to approach the mass and 
the temperatures of the non-relativistic limit to be dif- 
ferent (see below for more details). 

B. Non-relativistic behaviour 

As soon as the temperatures of the fermion fluids ap- 
proach the rest mass energy the scaling of the temper- 
atures and chemical potentials will start to be different. 
In the following equations n, T and /i refer to particles or 
antiparticles, not both. From eqs. (fT3|) and (f26|) we see 
that in the non-degenerate limit 

.«C,T3/^exp(^)=^ (30) 

If we now assume that ^ — m and T just scale like powers 
of the scale parameter 

fi-m = C2 Iti" and T = To ja^ (31) 

where Ci, C2, To, no, x and y are constants, we arrive at 

^ = Ciro'/V3-/2expga^-) (32) 

or 

a-^(ln^ + ^(,-2)lna).| (33) 

This can only be fuUfilled for arbitrary scales if a; = y = 2, 
meaning that 

T oc oT^ and (/x — m) oc oT^ (34) 

Note that this infers that the Boltzmann factor 
exp ((/i — m)/r) stays constant in the non-relativistic 
limit. The energy density and pressure will then scale 
as e « mn oc and p ~ Be'l^ oc in the non- 
relativistic degenerate limit. We stress that this is ex- 
actly the same scaling as that of a Boltzmann gas, where 
p « nT oc as T oc a~^. It is noteworthy that the 
evolution of energy density and pressure with the scale 
parameter does not change in both the ultrarelativistic 
and the non-relativistic limit when a finite chemical po- 
tential is added. 



C. Big bang nucleosynthesis 

The most important constaint on a chemical potential 
for dark matter can be deduced from big bang nucleo- 
synthesis (BEN) calculations. The abundances of light 



elements are sensitive to the expansion rate of the uni- 
verse during BEN. The expansion rate on the other hand 
is determined by the energy density of relativistic parti- 
cles (radiation domination) during this short epoch (a 
few MeV > Tphoton > 30 kcV) which is usually expressed 
in terms of the effective number of neutrino families N,y, 
for a detailed description see Ref. For standard cos- 
mology Ni, = 3 (or to be more precise = 3.046, which 
is due to the small overlap of neutrino decoupling and 
e~^e~- annihilation) while the values of deduced from 
comparison of the observed abundances of light elements 
to BEN-calculations and from WMAP data differ and 
have relatively large error bars. For example the au- 
thors of Ref. [iil find N,, = 3.08^};^^ at 2a error us- 
ing a combined ^He-abundance and CME results for the 
baryon abundance. On the other hand Ref. ^27,] arrives 
at = 5.2^2 2 '^'^ error using CME and large scale 
structure data. Ref. '28*1 find a similar upper limit for 
= 3.2^2.1 at 2a error using WMAP and BEN includ- 
ing ^He- and deuterium-abundances. In a more recent 
paper de Bernardis et al. [29^ find = 3.7^11 at 2a 
using CMB data, LSS, Supernova data and other inde- 
pendent measurements of the Hubble parameter, which 
is more restictive than Ref. (27l . [28j . So although N^, is in 
any case compatible with the standard value of N^, — 3 
it can be stressed that some of the observations seem to 
favour a larger value, as for example found in Ref. (27j by 
including large scale structure data. We will give upper 
and lower bounds on the chemical potential and particle 
mass to allow comparison to these different results. The 
standalone BEN results from Cyburt et al. [2^ are nev- 
ertheless the most relevant (and most restrictive) for our 
approach since all the other measurements of N^, refer 
to later stages in the evolution of the universe at which 
our dark matter particle will already be non-relativistic 
and thus does not contribute to the effective number of 
relativistic particle species. 

Next we will estimate the allowed window for the value of 
the chemical potential of the dark matter particle. One 
has to be careful with the definition of at this point 
since the important quantity is the temperature of the 
neutrino distribution, which will differ from the photon 
temperature at BEN. This is due to the annihilation of e+ 
and e~ below 0.5 MeV, which occurs after the decoupling 
of neutrinos at around 1 MeV, and the subsequent release 
of entropy almost completely into the photon fluid. So 
we flnd the following constraint on the energy density of 
a fermionic dark matter particle Ae^ei > f^DM- 

^^"120^'' - ■ 1,120^^^^ 



where T^ is the neutrino temperature. Using eq. ^ and 
deflning 
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(36) 



106.75 



2+ I (2 -2 + 2 -3) 



4/3 



21.3 



for decoupling when all standard model particles are 
present and in equilibrium. We arrive at the following 
constraint on fioM /Tdm- 

9dm ITT'' \Tdm J iTT^ \Tdm J 

Solving for the largest value of udm/Tdm that still full- 
fills this relation one arives at 



fJ'DM 



Tdm 



\ 




(38) 



This constraint is of course only valid for fermions that 
are still relativistic at BBN. It does not apply for particles 
that were non-relativistic at decoupling or non-thermally 
produced. This is due to the fact that the contribution 
of non-relativistic particles to the total energy density at 
BBN must be negligible not to overdose the universe to- 
day, consequently one can hardly probe udm /Tdm for 
such particles via BBN. 
The constraints on the maximum allowed values for 



to 2). For AA^^ = 0.094 a non-zero chemical potential 
can be excluded, for AA^;^ < 0.094 even a non-degenerate 
fermion family can be excluded (the T^^.^-term in eq. ([35|l 
always remains!). This limit will of course be shifted to a 
higher value of AiV^ if g^j^ is smaller at decoupling than 
our reference value of 106.75. For AiV^"°^ = 5 0, [111 
the maximum allowed value for the chemical potential is 
6.3 times the dark matter temperature. Futhermore even 
for AiVy — 1.63 [26], the ratio udm /Tdm can be as large 
as 4.4. Figure[l]also shows the allowed chemical potential 
for J = 50 at DM-decoupling and for hypothetical val- 
ues of 200 and 1000 showing that even such an enormous 
increase in g'^^^j will allow a maximum chemical potential 
of only a factor of about 2 larger compared to decoupling 
while gljj = 106.75. 

In principle g^^j j{TDec) could be much higher than the 
standard model value for example due to numerous addi- 
tional super symmetric particles. E.g. in a gas of strings 
in the very early universe there are exponentially many 
states present (see for example Ref. [30, 31]). As we will 
see later very large values of glfj{TDec) are essential for 
this approach to work. 

Of course there are other mechanisms that could lead 
to a variation in the extracted effective number of neu- 
trino families from the standard value of = 3.046. 
For example an evolving gravitational constant G could 
also lead to a difference in the expansion rate at BBN 
(see e.g. Ref. [11]), as it could also be the case in cosmo- 
logical models with more than the standard 4 spacetime 
dimensions (see e.g. Ref. jssj). 
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FIG. 1: Maximum allowed dark matter chemical potential 
over temperature as a function of the number of additional 
effective neutrino families for different values of the effective 
number of degrees of freedom gl^f at decoupling. The differ- 
ent lines correspond to the cases glf fiToec) ~ 1000 (dashed), 
200 (sohd gray), 106.75 (solid) and 50 (dash dotted) 

fJ-DM /Tdm are plotted vs. AN^, in figure [T] (g/^M is set 



V. NUMERICAL RESULTS 



In the following we examine one specific parameter set 
of a flat universe {^tot = 1) in our numerical treatment 
with the present day values = 0.72 (simple vacuum 
energy), Q^^i = 0.236, = 0.044 and n°^^ = 8.5 ■ 
10~^, similar to the parameters obtained from WMAP- 
data using the ACDM-model [IJ. ilx as usual denotes 
the fraction of the critical energy density of the respec- 
tive energy component today. We use a Hubble param- 
eter of the value Hq = 70 km s^^ Mpc^^. We start our 
calculation of the Fermi-integral-sets (e,p, n) and (e,p, n) 
from eqs. (|4]), (O, ([6]) when the particles and antiparti- 
cles are still ultra-relativistic and use number and energy- 
momentum-conservation, eq. (jl7|) , for both fermions and 
antifermions separately to find the scaling of the chem- 
ical potentials and temperatures. When the latter have 
reached non-relativistic scaling behaviour with sufficient 
precision we continue using the non-relativistic expres- 
sion of eq. p^ . 

We add photons, three massless non-degenerate neutrino 
families and baryons as particle background. To find the 
time dependence of the scale parameter a we use the 
first Friedmann-equation for a flat universe (i.e. with zero 
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curvature): 



SttG 



^tot 



(39) 



goM is in this case equal to four since particles and an- 
tiparticles contribute equally to the energy density. The 
lower limit on the mass (in other words AiV^ = AiV™°^) 
can be estimated by: 



Finally we calculate the free streaming length for 
fermions and anti-fermions to find additonal constraints 
on the particle mass and show that the scaling of the 
Jeans length and Jeans mass is exactly that of ordinary 
warm dark matter (WDM). 



A. Mass range 



3-Tl 



,0 / Tdm £7+1/ 



1.6 



"dm ^2 



106.75 



glffiToec) 



7 ^UM / TDM>m '""l+v 
1/3 

( 9DM \ 



AN, 



"dm 



0.236 / V0.70 



V 2 / 

2 

eV 



For a given initial value of fi /Tdm wc then use the 
condition ri^JAf ^ 0.236 to pin down the particle mass 
which is sufficient to close the universe. Figure [5] shows 
the allowed range of the mass, assuming AiV™°^ = 5 (i.e 
f^DAi/TDM < 6.3), were we find 1.8 eV < m < 53 eV. For 
^j\^max = X.63 (or a maximum initial value of ^dm /Tdm 
— 4.38) the corresponding mass constraint is 4.4 eV < 
m < 53 eV. For small initial values of (J-dm /Tdm one 
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FIG. 2: particle mass sufficient to close the universe for the al- 
lowed range of the initial chemical potential over temperature 
ratio 

finds the expected T'^-plateau due to the first term in 
the ultra-relativistic energy density, see eq. ([H), which is 
independent of the chemical potential. The critical mass 
is determined just by the relativistic degrees of freedom 
at decoupling and the expected density today (or at any 
other arbitrary point in the non-relativistic regime). A 
good estimate of the upper limit on the mass (meaning 
for udm/Tdm ~ 0) is 
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where in both cases we assume that the fermions be- 
come non-relativistic at T ss m/3. Most of the parame- 
ter range is dominated by the contribution of the second 
term in eq. ([5]), only for the values above udm/Tdm ~ 
4.4 the third term oc /z^ becomes dominant. For a larger 
effective number of degrees of freedom at decoupling the 
mass window is shifted to higher masses by a factor of 
glf f{TDec) / glf f{SM), where gljf{SM) corresponds to 
all standard model particles. In this case the highest al- 
lowed value of the initial chemical potential will be larger 
(see section IIV C[) meaning that the size of the allowed 
mass window will grow. It is also interesting to notice 
that the lower limit on the mass only has a slight de- 
pendence on the degrees of freedom at decoupling and is 
therefore always close to IcV for all discussed values of 

glff{TDec)- 

In figure [3] we show the dependence of the particle mass 
on the degrees of freedom at decoupling for initial val- 
ues of n/T = 0.1,1.0,2.0, which is the most interesting 
range of parameters as we shall see later when discussing 
structure formation. 
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FIG. 3: darli matter particle mass sufficient to close the uni- 
verse vs. the degrees of freedom at decoupling for several 
values of the initial chemical potential over temperature ratio 
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B. Scaling of the temperature and chemical 
potential 

Figures [4] and [5] show the scahng of the dark matter 
temperature and chemical potential for initial values of 
IJ-dm/Tdm = 5.5, 0.1 and a fermion mass of 2.54 eV and 
53.3 eV, respectively. In the ultra-relativistic regime (i.e. 
for small values of a) we recover the simple 1/a - scaling 
for /i DM and T^m as expected from eq. (1281) . In the non- 
relativistic regime we find that both chemical potentials 
are aproaching the mass and the temperatures scale as 
a^^ like in eq. (|13p . Furthermore, one can see that the 



tory for number conservation of both species (in figure [5] 
the difference is too small to be visible) . The chemical po- 
tential of the antifermions, /2, approaches the mass from 
negative values with a similar rate as /i making it look 
like it is shooting up exponentially in the log- log-plot. 
In figure [5l the case for small initial fJ-DM/ToM, we see 
that both the particle and the antiparticle chemical po- 
tential are temporarily below the mass and change their 
scaling only when the temperature approaches the mass. 



C. Equation of state 
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T antifermion 
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FIG. 4: Temperatures and chemical potentials for an initial 
fi/T = 5.5 and a mass of 2.54 eV 
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5: Temperatures and chemical potentials for an initial 
= 0.1 and a mass of 53.3 eV 



temperature of fermions will be higher than that of an- 
tifermions in the non-relativistic regime, which is manda- 



As outlined in sections IIII Al and IIII Bl one expects a 
radiation-like EoS in the ultra-relativistic limit as well 
as a polytropic equation of state with an exponent of 
5/3 in the non-relativistic limit for both a degenerate 
and a non-degenerate Fermi gas. Actually one does not 
expect to see the non-relativistic degenerate limit since 
the highest allowed initial value of IJ^dm /Tdm of 6.3 is 
only slightly higher than the value at which the //''-Term 
in ([8]) becomes larger than the /i^T^-Term at a value of 
4.4. So we actually examine only the medium- to non- 
degenerate region of a free Fermi gas. In Figure [S] the 
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FIG. 6: EoS for an initial /i/T = 5.5 and a mass of 2.54 eV 

EoS of fermions and antifermions for an initial /i/T = 5.5 
and a mass of 2.54 eV is shown. The EoS differs in the 
non-relativistic regime for a non-zero initial [Ibm ITbm 
because the fermions are (at most) semi-degenerate 
while the antifermions are non-degenerate and very di- 
lute (in comparison to the fermions). When comparing 
the fermion-pressure and the antifermion-pressure in the 
non-reltivistic regime at the same energy density for each 
particle-distribution the antifermions have a higher pres- 
sure because their temperature is higher for the same 
value of the energy density. The pressure of fermions is 
of course always larger than the pressure of antifermions 
at the same scale parameter for positive initial fermion 
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chemical potential, because the energy density of an- 
tifermions can in fact be orders of magnitude smaller. 



D. Scaling of energy densities, number densities 
and pressures 

For a flat universe the total density is equal to the crit- 
ical density at all values of the scale parameter etot (a) — 
^crit (a) ■ This can be used to deduce the scaling of the in- 
dividual normalized energy densities (usually denoted by 
fi). For example the normalized radiation energy density 
is given by 



£7+1/ 



(40) 



where e is the energy density of fermions and e is the 
one of antifermions. We show the scaling of the mean 
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FIG. 8: As fig.Elbut for an initial fi/T = 0.1 
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FIG. 7: Fractions of the critical density for an initial fi/T — 
5.5, photons and neutrinos (solid gray), dark matter(fermions: 
solid, antifermions: dash dotted gray), baryons (dash dotted) 
and dark energy (dashed) 

energy densities normalized to the critical density in fig- 
ures [7] and [8] for the cases udm/Tdm — 0.1 and 5.5. 
The plots show nicely what is often called the "coinci- 
dence problem" , namely that dark energy seems to to 
become dominant just below a redshift z = ^ — 1 « 0.5. 
On can see that the contribution from fermions and an- 
tifermions is constant during radiation domination until 
they become non-relativistic and can constitute a sub- 
stantial fraction of the total energy budget in the very 
early universe. Then their normalized energy densities 
grow linearly with a until non-relativistic particles domi- 
nate. In contrast to that cold dark matter would in these 
plots scale just like the baryon density and exactly coin- 
cide with the scaling behaviour of the sum of the fermions 
and antifermion contributions in the matter dominated 
regime, i.e. the CDM contribution is highly suppressed 
at early times. 



In Figure [9] we show an example of the scaling of en- 
ergy density, pressure and number density of fermions 
and antifermions an initial (J-dai /Tdm = 5.5. In the 
ultra-relativistic limit we recover radiation like scaling of 
energy density and pressure being proportional to a"**, 
in the non relativistic limit the energy densities scale 
like matter (e cx a^^) and the pressure drops rapidly 
{p cx a~^). The number densities always have a slope of 
—3 as demanded by number conservation of particles and 
antiparticles. 
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FIG. 9: Evolution of the energy density (solid), pressure 
(dashed) and number density (dash dotted) of fermions 
(black) and antifermions (gray) for an initial fi/T — 5.5 and 
a mass of 2.54 eV 



E. Age of the universe 

It is generic of the cases studied (r2° = 0.72, = 
0.236) to produce an almost fixed age of the universe. 



10 



independent of the initial chemical potential over tem- 
perature ratio of the dark matter component. This is 
due to the fact that the dark matter has to become non- 
relativistic before the decoupling of photons and baryons 
occurs at around 380 000 years after the big bang. So 
any increase in the energy density at an earlier stage will 
only reduce the age of the universe by less than 380 000 
years due to an increased Hubble parameter at earlier 
times. 

So for nl = 0.72, n%j^i = 0.236 we find the age of the 
universe is the well known 13.7 Gyrs (as in ACDM), for 
example in a flat Sl^/ = 1 universe the age would be only 
9.3 Gyrs(!). The latter case (although already being ex- 
cluded by WMAP, if the spectrum of initial fluctuations 
is a simple power law fT] ) has the common problem of all 
flat matter dominated models to be in violation of some 
independent lower limits on the age of the universe. This 
is due to the today well messured value of the bubble pa- 
rameter of around 70 km s~^ Mpc^^. For example from 
stellar evolution [sj] one finds a lower limit of 10 Gyrs on 
the age of uj Gen and Ref. [s^l finds M92 to be at least 13 
Gyrs old. From the cooling of white dwarfs in globular 
clusters [s^ derives an age of more than 10.3 Gyrs (even 
at a 2a lower limit) on the age of M4, while in [331 NGG 
6397 is found to be 11.47±0.47 Gyrs. 
A significantly lower value of the bubble parameter [h < 
0.5) would be needed to solve this age problem and reach 
an age of more than 13 Gyrs for a flat matter dominaded 
universe also. 



Structure formation 



which is normaly in between the scale where particles be- 
come non-relativistic and the scale where matter starts 
to dominate. The mean velocity of fermions can be cal- 
culated from the distributions function 



(v) 



goM 



1 



E{k) exp(^^) + 1 



(43) 



The corresponding mass enclosed in a sphere of radius 
\fs is the free streaming mass: 



47r 3 



(44) 



From the comparison of numerical n-body simulations of 
warm dark matter and cold dark matter Ref. (40| claim 
that a (thermal) warm dark matter particle of a mass be- 
low 750 cV would be incompatible with observations of 
the Ly a forest at a redshift of z « 3. On the other hand, 
Boehm et al. find with higher resolution simulations 
that warm dark matter particle masses down to 600 eV 
would produce a power spectrum that is virtually indis- 
tinguishable from a simple cold dark matter spectrum be- 
low z ~ 2. This effect arises from the fact, that as soon as 
large scale fluctuations go non-linear, small scale fluctua- 
tions will start to grow non-linear as well and quickly re- 
generate an initially suppressed power spectrum on small 
scales. The conclusion of [4l| is that the primordial power 
spectrum could have been exponentially surpressed ini- 
tially up to mass scales of ~ IO^Mq without being ob- 
servable in large scale structure surveys today. We will 
now use this upper limit on the free streaming mass to 
find a lower limit on the particle mass or the degrees of 
freedom at decoupling using the prescription of eq. (|42p . 



Now we will discuss more restrictive limits on the parti- 
cle mass and degrees of freedom at decoupling which can 
be deduced from the shape of the matter power spec- 
trum. In the second part of this section we will show 
that the general scaling of the relevant quantities for lin- 
ear structure formation do not change when adding a 
finite chemical potential. 



1. Free streaming 

For coUisionless particles (that have decoupled while 
being ultra relativistic) the most important damping 
scale is the free streaming length, as found by (sl]. The 
comoving free streaming length is given by 



A 



FS 



MO) 

a{t') 



dt' (assuming tjjec t) (41) 



which is basically the distance a particle can travel in 
an expanding background since decoupling. Interestingly 
Bond and Szalay [39| found that the relevant comoving 
damping scale from free streaming of coUisionless parti- 
cles can be estimated by 



A_FS ~ 3.2 



{v{t)) t 



a{t) 



(42) 
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FIG. 10: Free streaming length of fermions vs. the degrees of 
freedom at decoupling for several values of the initial chemical 
potential over temperature ratio 

In Figure 10 we show the resulting free streaming 
length and free streaming mass of fermions as a func- 
tion of the degrees of freedom at decoupling for several 
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FIG. 11: Free streaming mass of fermions vs. the degrees of 
freedom at decoupling for several values of the initial chemical 
potential over temperature ratio 
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FIG. 12: Free streaming mass vs. the fermion mass for several 
values of the initial chemical potential over temperature ratio 



magnitude. [As a side remark: of course for a Majorana 
fermion glfj{TDec) does only need to be half as large to 
reach the same mass as for a Dirac fermion with /i/T k, 0] 
So structure formation gives far more restrictive 
bounds on the initial /i/T, namely that the ratio should 
not exceed unity by far (keeping in mind that BBN al- 
lowed /i/T to be up to 14 for glffiToec) = 10^). Nev- 
ertheless one can in principle only give a lower bound 
on g^j:j:(Toec) since the free streaming mass is mostly 
fixed by the particle mass (as one can see in Figure 
I12p . which increases linearly with .g|^j(rDec)j meaning 
that higher values of fi/T could be possible for even 
higher values of the degrees of freedom at decoupling 
{glffiTnec) > 4000). The sHght decrease in the free 
streaming mass for rising ^ldm /Tdm at fixed mjjM is 
caused by the increase in g'^j:j:{T]jec) in order to keep 
^'dm constant. A higher ^^^^(T^iec) means a lower tem- 
perature relative to the temperature of equilibrium par- 
ticles, meaning the dark matter particles become non- 
relativistic at a smaller scale parameter unr- Remem- 
bering eqn. ^ this means 



aNR oc 



T / TDM>m 



Of glffiToec) 



-1/3 



(45) 



at a fixed dark matter particle mass, which explains the 
small reduction in Mps- 

If the dark matter distribution should not exactly coin- 
cide with the baryon distribution up to such small scales 
as measured by the Ly a data (possibly measurable by 
future weak lensing surveys) then this could of course 
relax the given restrictions. 

Another way of giving a lower limit on the mass of 
dark matter candidates is the dark matter phase space 
density in galaxies. A more realistic version of the origi- 
nal Tremaine-Gunn-limit [4^ takes into account changes 
of the initial particle distibutions functions during col- 
lapse to a galaxy halo as discussed by Madsen [4^]. He 
finds a lower limit on the mass of a fermionic dark mat- 
ter particle of ^ 380 eV, which is less restrictive than the 
limits from the shape of the matter power spectrum. 



values of the initial ratio of ^i/T. The free streaming 
length and mass is the same for antifermions at vanish- 
ing chemical potential and slightly lower for increasing 
chemical potential because antifermions will then become 
non-relativistic earlier as we have seen in Figure 4. Since 
the relative number of antifermions is also suppressed in 
the latter case the total damping scale should still be 
well described by the damping scale of fermions. One 
can also see that even for a n/T = 2.0 one can stay 
below the threshold of ~ IO^Mq although this requires 
more than 2000 degrees of freedom at decoupling. Inter- 
estingly, even for a vanishing chemical potential a value 
of gl^jiToec) ~ 10^ (corresponding ttidm ~ 500 eV sim- 
ilar to the lowest mass value used in Ref. [4l|) is needed 
to be in agreement with Ly a measurements, which al- 
ready exceeds the standard model value by an order of 



2. Linear structure formation 

The simplest approach to describe the onset of struc- 
ture formation after matter radiation equality is to look 
at the linear equation for the density contrast of mat- 
ter 5{t, x) in comoving coordinates {x — r/a) which is 
defined by p{t,x) — p{t) {I + S{t,x)). p is the mean 
mass density of the corresponding matter component. In 
fourier space the evolution of the density contrast with 
time can be found by solving 

5k + 2H5k + - ^irGp j dk = (46) 

where 5k{t) corresponds to the density contrast for one 
specific wave number k describing density pertubations 
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at the comoving lengthscale Xcom = 27r/fc or the physical 
lengthscale Xphys — 2T:a/k. For a non-rclativistic col- 
hsionless fluid the velocity dispersion (w^) replaces the 
square of the speed of sound of a coUisional fluid and 
is in our approach given by (see Ref. [44] and section Hn]) 

(47) 

^ ' mn 5 op \teqj 
On the other hand 



for a coUisional non-relativistic fermion fluid. One can see 
here that the velocity dispersion (or the speed of sound) 
drops like for an ordinary warm dark matter scenario 
(Ref. [SElilillil) contrary to the generalized Chap- 
lygin gas which has a late time growing speed of sound. 
Now inserting the velocity dispersion into eq. (|46p one 
finds the critical physical wavelength, the Jeans length, 

The corresponding Jeans mass reads: 

(49) 

Hence, the scalings derived here turn out to be the same 
as for standard WDM or CDM. In principle one has 
to do a full calculation of the linear power spectrum of 
density fluctuations in this model, but we expect results 
similar to warm dark matter models, although the final 
shape of the spectrum may be a bit more complicated 
due to the interplay of two (more or less) different dark 
matter components. 



VI. CONCLUSION 

We have shown that a relic fermionic particle family 
with sizeable chemical potential cannot be excluded even 
from the most recent BBN-calculations and large scale 
structure data. Using different constraints on the energy 
density at BBN we find that even for decoupling when 
glff{TDec) is equal to the maximum standard model 
value of 106.75 the initial ratio of ji/T could have been as 
large as 4.4 (when using the results of ]26|]), 4.5 (for the 
limits of [1^) or 6.3 (according to [27l l28j|l. In this case 
the resulting mass range is at most 1.8 eV < m < 53 eV, 
where the upper bound scales linearly with glff{Tuec)- 
In figure [Ql we show the combined results on gljj{Tnec) 
and iriDM from floM, the free streaming mass Mps 
and the Tremaine-Gunn-limit up to very high values of 
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FIG. 13: Qef fiToec) vs. niDM parameter plane. Shaded re- 
gions are excluded by present day energy density (diagonal 
pattern), by the Tremaine-Gunn-limit (horizontal pattern) or 
by a too high free streaming mass (vertical pattern). 



glffiToec) — 5 • 10"*. The region with the diagonal pat- 
tern is excluded because of a too high present day energy 
density {^dm > 0.236 for h = 0.7) even for an initial 
/i/T = 0. The border line is just the case of vanishing 
chemical potential (compare to eqn. ((40)) ') and therefore 
also applies to an ordinary warm dark matter relic. For 
all points below this line there is always one combina- 
tion of muM, glffiToec) and n/T that gives the right 
present day energy density. The dashed lines show these 
right combinations for several initial values of /x/T from 
1 to 6. To make this clear, a (ttidm, glj fiToec)) point on 
the fi/T = 1 line gives the right ^dm — 0.236 today for a 
fermion with /i/T = 1 and a too high Q.dm for p,/T > 1 
and vice versa for ji/T < 1. A free streaming mass of 
more than ~ 10^ solar masses is excluded for warm dark 
matter as found by [4l[ . The region excluded by this con- 
straint is the area with the vertical pattern. Note that 
the limiting line has been extrapolated beyond the U, dm 
limit and values of /i/T above 6. The region with the 
horizontal pattern is excluded by the Tremaine-Gunn- 
limit on the phase s pac e density of fermions in galaxies 
as found by Madsen [4^ , which is independent of the de- 
grees of freedom at decoupling. Fermion masses of less 
than 380 eV are excluded, but as one can see this limit 
becomes relevant only for (?|yy(TDec) > 10^. The con- 
straints from BBN are actually so much less restrictive 
that they lie way outside of this plot. Remember that 
for g^j:j:{TDec) = 106.75 the upper bound on the initial 
p/T was 4.4 (for AiV^ < 1.63) and for glff{TDec) = 10^ 
it would be even p/T < 10.7. Now note that the limit on 
the free streaming mass rules out a non-zero p/T even for 
glffiToec) = 10'^. This case would require for example 
only a tiny additional ANi, « 5 • 10^"^ to the number of 
effective neutrinos families at BBN. The BBN constraints 
become even less restrictive with growing (?|f j-(Tbec) and 
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thus do not appear in this summary plot. As one can see 
from eqn. psp the maximum allowed value of ^l/T grows 
as {glfj{TDec)Y^^ for large values of glff- So to summa- 
rize these results, we only find absolute lower limits on 
the fermionic dark matter mass and the effective degrees 
of freedom at dark matter decoupling. The initial ratio 
of chemical potential over temperature can be bounded 
from above for a given value of g^ffiToec) but no abso- 
lute upper limit has been found. 

We have seen that at least an order of magnitude more 
effective degrees of freedom than the standard model 
value are needed to allow a fermionic dark matter parti- 
cle to be the dominating dark matter component (even 
with vanishing chemical potential) and this is also the 
case for an ordinary warm dark matter relic. In a gas of 
strings in the very early universe there are exponentially 
many states present (see for example Ref. [s^, [Ml), so 
this would be a possible scenario to account for the very 
high values oi g^j f[Tjjf,c) that are needed. 

It is noteworthy that the limits on m dm from the free 
streaming mass and from the phase space limit cross, 
implying that the two main problems of CDM that stan- 
dard warm dark matter can probably only resolve with 
two different WDM particle masses (i.e. the excessive 
amount of small structures and the cuspy core issue in 
standard CDM, see e.g. Ref. [4?|) may be resolvable in 
the fermionic dark matter approach without this discrep- 
ancy. 

Note that decoupling at g^j^iTjjec) > 100 needs a 
super-weak interaction between dark matter and stan- 
dard model particles, therefore constraints from cooling 
of supernovae and red giants that apply for light weakly 
interacting particles or axions do not apply here. 

We have calculated the scaling of temperature, chem- 
ical potential, equation of state, energy density, number 
density and pressure of fermions and antifermions and 
have shown that even their temperatures and equations 
of state (as well as all other mentioned quantities) will 
be different in the non-relativistic regime for an initially 
non-zero chemical potential. We demonstrated that the 



velocity dispersion, the Jeans length and the Jeans mass 
all have the same scaling as for regular cold or warm dark 
matter after adding a finite chemical potential. 

An interesting feature about the model studied is prob- 
ably that two different dark matter components can 
emerge from the same source, that do not necessarily 
clump with the same magnitude in the same region of 
space. This is due to a different temperature and veloc- 
ity dispersion in the non-relativistic regime. A detailed 
numerical simulation of structure formation in such an 
approach would be necessary to make a clear statement 
about differences in the matter power spectrum today 
compared to a warm dark matter particle with the same 
mass but zero initial chemical potential. 

Finally, it is important to keep in mind that it is of 
course assumed that the dark matter particle was in 
statistical equilibrium at some point in the evolution of 
the universe, otherwise it will not follow a Fermi-Dirac 
distribution. For example in the already mentioned 
string gas cosmology approach such an early period 
of equilibrium would naturally occur. No assumption 
was made regarding the production mechanism and the 
source of the matter-antimatter asymmetry. The scope 
of this work was primarily to determine if a fermionic 
dark matter particle could in principle have a sizable 
chemical potential using cosmological observations and 
not to motivate such a relic from an underlying theory. 
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